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Abstract 

We present a new method for modeling the hadronic recoil in W — > tv events produced at hadron colliders. 
The recoil is chosen from a library of recoils in Z — > ££ data events and overlaid on a simulated W — > iv 
event. Implementation of this method requires that the data recoil library describe the properties of the 
measured recoil as a function of the true, rather than the measured, transverse momentum of the boson. 
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We address this issue using a multidimensional Bayesian unfolding technique. We estimate the statistical 
and systematic uncertainties from this method for the W boson mass and width measurements assuming 1 
fb _1 of data from the Fermilab Tevatron. The uncertainties are found to be small and comparable to those 
of a more traditional parameterized recoil model. For the high precision measurements that will be possible 
with data from Run II of the Fermilab Tevatron and from the CERN LHC, the method presented in this 
paper may be advantageous, since it does not require an understanding of the measured recoil from first 
principles. 

Key words: W, Z, mass, width, hadron, collider, Tevatron, DO, recoil 
PACS: 12.15.Ji, 13.85.Qk, 14.70.Fm, 12.38.Be 



1. Introduction 

The W and Z bosons are massive gauge bosons that, along with the massless photon, mediate elec- 
troweak interactions. The predictions from the standard model (SM) of weak, electromagnetic, and strong 
interactions [l| for their masses and widths include radiative corrections from the top quark and the Higgs 
boson. When precision measurements of the W boson mass (Mw) are combined with measurements of the 
top quark mass and other electroweak observables, limits on the Higgs boson mass can be extracted. The 
W boson width (Tw) can be directly measured from the fraction of W bosons produced at high mass. It 
can also be inferred indirectly within the context of the SM from the leptonic branching fraction of the W 
boson. The branching fraction, in turn, can be inferred from the ratio of the W and Z boson cross-sections 
with additional theoretical inputs [2| . The direct measurement of Tw is sensitive to vertex corrections from 
physics beyond the SM. The current world average for Mw is 80.398 ± 0.025 GeV [3j and the current world 
average for Tw is 2. 106 ±0.050 GeV from direct measurements The large number of W bosons currently 
available in data samples collected at the Fermilab Tevatron collider and that will soon be available from 
the CERN LHC collider allow measurements of Mw and Tw with unprecedented precision provided the 
response of the detector can be modeled with sufficient accuracy. 

In pp and pp collisions, W and Z bosons are produced predominantly through quark-antiquark anni- 
hilation. Higher order processes can include radiated gluons or quarks that recoil against the boson and 
introduce non-zero boson transverse momentum Figure [T] shows an example diagram for the production 
of a W or Z boson with two radiated gluons in a pp collision. 
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We identify W and Z bosons through their leptonic decays (W — > iv and Z — > with I = e,fi) since 
these signatures have low backgrounds. The charged leptons can be detected by the calorimeter or the muon 
system, while the neutrino escapes undetected. We do not reconstruct particles whose momentum vectors 
are nearly parallel to the beam direction, and therefore we only use kinematic variables in the transverse 
plane that is perpendicular to the beam direction. The neutrino transverse momentum vector (pjf) is 
inferred from the missing transverse energy (ffi T ), which can be calculated using the transverse momenta of 
the charged lepton (p^) and the recoil system (ut)- 

$ T = -[pS + u T ]. (1) 

We measure ut by summing the observed transverse energy vectorially over all calorimeter cells that are 
not associated with the reconstructed charged lepton. 

The recoil system is difficult to model from first principles; unlike the decay lepton, it is a complicated 
quantity involving many particles, as well as effects related to accelerator and detector operation. The 
recoil system is a mixture of the "hard" recoil that balances the boson transverse momentum and "soft" 
contributions, such as particles produced by the spectator quarks (the "underlying event"), other pp (pp) 
collisions in the same bunch crossing, electronics noise, and residual energy in the detector from previous 
bunch crossings ( "pileup" ) . Figure [2]shows transverse energies recorded in the calorimeter of the DO detector 
versus azimuthal angle and pseudorapidity [f| for a typical W — > ev candidate. The diffuse energy deposits 
spread over much of the detector are due to the recoil system. 

The various components of this measured recoil system have different dependences on instantaneous 
luminosity. For example, pileup and additional inelastic collisions scale with instantaneous luminosity, while 
the contribution from the underlying event is luminosity independent. Moreover, detector effects such 
as suppression of calorimeter cells with low energy to minimize the event size (zero-suppression cuts) can 
introduce correlations between the calorimeter response to the hard component and various soft components 
in the event, so that the detector responses to these components cannot be modeled independently. 

Two approaches have been used previously to model the W boson event, including the recoil system. 

One method takes the underlying physics from a standard Monte Carlo (MC) event generator and smears 

it parametrically to reproduce detector effects . The parameters are tuned to an independent but 

kinematically similar data set, namely Z — » 11 events. The second approach, or "Ratio Method", constructs 

Mt template distributions by directly taking Z — ► It events from collider data, treating one of the leptons 

as a neutrino [10]. The ratio of the Z boson mass to the corresponding W boson mass, taken together with 

the precisely measured Z boson mass |ll| from the CERN LEP collider determines the W boson mass. In 

this method, small differences in the Z and W boson line shapes and transverse momentum and rapidity 

distributions of the decay leptons must be taken into account. 

This paper presents a novel approach for modeling the recoil system in W — > iv events at hadron 
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colliders that uses recoils extracted directly from Z — > 11 collider data. The Z — > ££ data provide a mapping 
between the Z boson transverse momentum (pj?) and the transverse momentum of the recoil system (ut)- 



Versions of the recoil library approach have been proposed in the past [12| that used the map between the 
reconstructed pj? and the measured ut directly. In this paper we use a two-dimensional Bayesian unfolding 
method to derive a relation between the true and the measured ut, which allows the simulation of the 
recoil system for the same generator level value of the W boson transverse momentum (py ). 

The recoil library method presented in this paper has many advantages. Since the recoils are taken 
directly from Z — > 11 data, they reflect the event-by-event response and resolution of the detector. The 
additional soft recoil is built in, as is the complicated zero-suppression-induced correlations between it and 
the hard component of the recoil. Proper scaling of the recoil system with instantaneous luminosity is 
automatic since the W and Z samples have similar instantaneous luminosity profile. The most significant 
advantage of this method lies in its simplicity since it does not require a first-principles understanding of the 
recoil system and has no adjustable parameters. The dominant systematic uncertainties of this approach 
come from the limited statistical power of the Z — > £1 recoil library, as is true with the other methods. 

In this paper, we outline the implementation of this method. The method is tested using the W boson 
mass and width measurements. Only the electron decay channel is discussed, but our method can also be 
used in the muon decay channel. The detector and selection criteria are described in Section [2] The MC 
simulation samples used are described in Section[3] We discuss the method in Section^ In Sections [5] and 
we assess the uncertainty on the W boson mass and width measurements, and compare the performance of 
this method with that of a parameterized recoil model. The paper concludes in Section [71 

2. The W and Z Boson Measurements 

We evaluate the recoil library method by estimating biases and statistical and systematic uncertainties on 
the W boson mass and width measurements in the electron channel. The test is performed using simulations 
of the Run II DO detector at the Fermilab Tevatron, a pp collider with center of mass energy s/s =1.96 TeV. 
Statistical uncertainties are estimated for a data sample corresponding to 1 fb _1 . 

2.1. The DO Detector 

The DO detector Q,Q] consists of a magnetic central-tracking system, comprised of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT), both located within a 2 T superconducting solenoidal 
magnet. The SMT and CFT cover \q D \ < 3.0 6] and \r) D \ < 1.8, respectively. 

Three uranium/liquid argon calorimeters measure particle energies. The central calorimeter (CC) covers 
\rio\ < 1.1, and two end calorimeters (EC) extend coverage to \rjrj\ ~ 4.2. In addition to the preshower 
detectors, scintillators between the CC and EC cryostats provide sampling of developing showers at 1.1 < 



\t)d\ < 1-4. The CC is segmented in depth into eight layers. The first four layers are used primarily to measure 
the energies of photons and electrons and are collectively called the electromagnetic (EM) calorimeter. The 
remaining four layers, along with the first four, are used to measure the energies of hadrons. Most layers 
are segmented into 0.1 x 0.1 regions in (77, <fi) [6] space. The third layer of the EM calorimeter is segmented 
into 0.05 x 0.05 regions. 

A muon system is located beyond the calorimetry and consists of a layer of tracking detectors and 
scintillator trigger counters before 1.8 T iron toriods, followed by two similar layers after the toroids. Tracking 
at \t)d\ < 1 relies on 10 cm wide drift tubes, while 1 cm mini-drift tubes are used at 1 < \rjn\ < 2. 

Scintillation counters covering 2.7 < \t)d\ < 4.4 are used to measure luminosity and to indicate the 
presence of an inelastic collision in beams crossing. 

2.2. Measurement strategies for jV% and Tw 

The W boson mass is measured from distributions of the following observables: the electron transverse 
momentum jjrf, the missing transverse energy tfJ T , and the transverse mass Mr, given by 



where A(f> is the opening angle between jjrf and Tp T in the transverse plane. The data distributions are 
compared with probability density functions from MC simulations generated with various input W boson 
mass values ("templates"). A binned negative log- likelihood method is used to extract Mw- The W 
boson width is measured using a similar method, except that only events in the high tail region of the Mj- 
distribution are used. For the mass measurement, the fit ranges we used are [30, 48] GeV for the \p^\ and 
$5 T | distributions, and [60, 90] GeV for the Mr distribution. For the width measurement, we fit the Mr 
distribution over the range [100, 200] GeV. 

2.3. Selection criteria 

A W boson candidate is identified as an isolated electromagnetic cluster accompanied by large \$ T \. The 
electron candidate is required to have a shower shape consistent with that of an electron, \p^\ > 25 GeV, 
and |?7r>| < 1.05. To further reduce backgrounds, the electron candidate is required to be spatially matched 
to a reconstructed track in the central tracking system. Additionally, we require |^ T | > 25 GeV, \ut\ < 15 
GeV, and 50 < Mr < 200 GeV. Z boson candidates are identified as events containing two such electrons 
with di-electron invariant mass 70 < M ee < 110 GeV and |-ur| < 15 GeV. The selection on \ut\ helps 
to suppress background and to reduce the sensitivity of the measurement to uncertainties on the detector 
model and the theoretical description of the p^ distribution. Since the Z sample has fewer events at high 
p^,, the detector and theoretical models are best constrained at low boson px- 

For this analysis, both electrons from the Z boson are required to be in the central region of the 
calorimeter because the unfolding requires well-understood detector resolutions. 




(2) 
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3. MC Simulation Samples 



In this paper we use three different MC simulations. Two of these are fast MC simulations and the third 
is a detailed full MC simulation using GEANT 15| . The two fast MC simulations are built around a common 



event generator and parametric model for the electron measurement, but with different recoil models. One 
uses a traditional parameterized method to model the recoil system, which we call "the parameterized recoil 
method" . The other uses our new method, which we call "the recoil library method" . The full MC Z — > ee 
sample has the equivalent of 6.0 fb -1 in integrated luminosity, and the full MC W — > ev sample corresponds 
to 2.5 ft) -1 . 

For both fast MCs, the pythia [l(J event generator is used to simulate the production and decay of the W 
boson, as well as any final state radiation (FSR) photons. FSR photons, if sufficiently close to the electron, 
are merged with the electron. After the event kinematics are generated at the four-vector level, detector 
efficiencies and energy response and resolution for the electron are applied. These parameterizations are 
measured using Z — > ee events from either collider data or full MC, depending on the study. A parametric 
energy dependent model for resolution effects is used. Parameterized efficiencies for data selection are 
prepared for comparing with either data or full MC as a function of electron [p<j?|, rj e , the component of 
the recoil along the electron direction, the total hadronic activity in the event, and the reconstructed z 
coordinate along the beam line where the hard scattering occurred. The recoil system is then modeled 
either using the recoil library or the parameterized model. 

The parameterized recoil method models the detector response to the hard recoil using a two-dimensional 
parameterization of the response (both magnitude and direction) as estimated using GEANT-simulated Z — > 
vv events. The underlying event is modeled using the measured Jfi T distribution from data taken with a 
trigger that requires energy in the luminosity monitors ( "minimum bias events" ) , and pileup and additional 
interactions are modeled using the measured Jfi T distributions from unsuppressed data taken on random 
beam crossings ("zero bias events"). These are combined with the hard recoil, and data-tuned corrections 
are applied to account on average for correlations between the "hard" and "soft" recoil. The correction 
parameters are tuned to Z —* ee control samples. The parametric methods of modeling the recoil are 



further discussed in Refs. [17( and resemble approaches used in earlier DO and CDF measurements at the 
Tevatron ^, 3, Q . The recoil library method of modeling the recoil is discussed in detail in Section [U 

The GEANT-based MC simulation also uses PYTHIA to simulate the production and decay of the W 
boson, as well as the underlying event and any FSR photons. These events are then propagated through 
a detailed description of the detector. Zero bias collider data collected by the DO detector with a similar 
instantaneous luminosity profile as the W — > ev collider data sample are overlaid on the full MC simulation 
to model additional collisions and noise in the detector. These events are processed through the same full 
set of DO reconstruction programs as data. 



4. The Recoil Library Method 



4-1. Overview 

The recoil library is built from Z — ► ee events. Because the electron energies and angles are well 
measured, the measured p T z from the two electrons provides a good first approximation of the true Vt ■ 
An unfolding procedure allows the transformation of the two-dimensional distribution of the measured |p T z | 
and measured \ut\ to that of the true \p^\ and measured \ut\- The opening angle between the measured 
PjF and the measured ut is also unfolded to the opening angle between the true p^ and the measured ut 
during this procedure. A map between the true \Pt\i the measured \ut\, and the scalar Et (SEt), which 
is defined as the scalar sum of the transverse energies of all calorimeter cells except those that belong to 
the reconstructed electrons, is also produced. This map is not used by the recoil model, but is needed by 
the electron efficiency model. The final result of the recoil library is the ut for an event, referenced to the 
true Pt . This vector substitutes for the equivalent vector obtained in the parametrized recoil model. All 
further corrections for efficiencies due to the recoil system are the same for both the recoil library and the 
parametrized recoil model. 

Figure [3] shows some examples of the distribution of the component of the measured recoil along the Z 
boson direction and perpendicular to the Z boson direction. 

4-2. Preparing the recoil library 

Before producing a binned recoil library, certain cvcnt-by-cvcnt corrections must be applied to the mea- 
sured recoil system. We need to remove any biases in the measured recoil distribution due to the Z boson 
selection requirements. Electron identification requirements, for example, preferentially reject events with 
significant hadronic activity. Events with significant hadronic activity also have poorer recoil resolution than 
events with little hadronic activity. Since Z bosons contain two electrons while W bosons only have one, the 
bias will not be the same. The electrons from Z boson decays also have a higher average \pjt\ and a different 
rf distribution than electrons from W boson decays. To account for this, we remove the biases from the Z 
boson selection, and then, when a W candidate is made using the recoil library, the biases appropriate for 
a W candidate are added, as described in Ref. In this section, we describe these corrections to the Z 
boson sample. 

4-2.1. Removing the two electrons from Z boson events 

The recoil system for Z — > ee events is defined as the energies in all calorimeter cells excluding those that 
belong to the two electrons. Since the recoil system will in general deposit energy in these cells, excluding 
them biases the component of the recoil along the electron's direction. We correct this effect by adding back 
an approximation of the underlying energy. 

9 



This correction (denoted by Aitjj) depends on uf, (the projection of ut along the electron transverse 
direction), instantaneous luminosity, and electron rj e , and is estimated using the energies deposited in equiv- 
alent cells that are separated in <j) from the electron in W — > ev events. In addition to correcting for the 
recoil energy under the electron cluster, we also correct for electron energy that leaks out of the cluster. For 
Z boson events, these corrections are made for both electrons. In Section [5] we estimate the uncertainty due 
to these corrections. 

4-2.2. Minimizing the effects of FSR photons 

The full MC simulation indicates that roughly 6% of the Z — > ee events contain FSR photons with 
E T > 400 MeV that are sufficiently far from the electrons that the electron clustering algorithm at DO 
does not merge them with a reconstructed electron. These photons are thus incorrectly included in the 
measurement of Ut, instead of in pj? , resulting in a correlated bias. Since Z — > ee events contain more 
FSR photons than W — > ev events do, the recoil library built using Z bosons will contain on average larger 
contributions from FSR photons. 

Ideally, these FSR photons could be removed from the recoil file, and the effect could be separately 
modeled within the fast MC simulation. Since it is difficult to identify these FSR photons on an event-by- 
event basis, the effect is reduced by raising the lower limit on the effective reconstructed di-electron invariant 
mass to 85 GeV, reducing the fraction of events with a high Et FSR photon by 25%. 

The effect of the remaining photons is small because, for a low px W boson, Mt ~ 2|p* 3 ?| + ujj. Therefore, 
the photons will create a bias on the mass only if they produce a bias in the component of ut parallel to 
the electron direction. While the overlaid recoil is rotated so that the direction of its corresponding Z boson 
matches that of the simulated W boson, the directions of the decay electrons from Z and W are largely 
uncorrelated, and the bias is mostly canceled for measurements using the Mt spectrum. In Section [5] the 
bias due to the FSR photons is estimated. 

4-2.3. Correcting for electron selection efficiencies 

The selection criteria for W and Z candidates can introduce biases between the electron and the recoil 
system. Since the kinematic and geometric properties of W candidates are not identical to those of Z 
candidates, they have different biases. 

The two components of the electron selection efficiency model that most strongly affect these biases are 
the SEt efficiency and the itjj efficiency. The SEt efficiency describes the electron identification probability 
as a function of the overall activity in the detector. The ujj efficiency describes the probability of electron 
identification as a function of ujj. This probability decreases with increasing hadronic activity along the 
electron direction. 

Since the recoil library is built from Z — » ee events, we need to remove the biases introduced by the 
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selection requirements on the two electrons. We correct for the efficiencies by weighting each event in the 
Z boson recoil library by l/ e "|j x ^/ e SE T , where e u jj is the product of the u| efficiencies and esE T is the 
product of the SEt efficiencies for the two electrons in each Z candidate. 

When W boson events are produced from a fast MC using the recoil library, the map between the true 
Pt , measured ut, and SEt is used to introduce the biases appropriate for W bosons from these efficiencies. 
To simulate a W boson event, a random recoil is chosen from the recoil library corresponding to the true 
W boson p^ , and a random SEt is chosen from the SEt distribution corresponding to the true W boson 
p^ and the chosen recoil ut- The ujj efficiency and SEt efficiency are then applied to the electron from W 
boson decays. 

4.3. Unfolding method 

After the recoils have been corrected as discussed above, the transformation from measured pf 2 and 
measured ut to true Pt and measured ut is done using a Bayesian unfolding technique. 

4-3.1. Multidimensional unfolding using Bayes's Theorem 

Unfolding is a mathematically challenging problem, since it involves the reversal of a random process. 
Because a given true state can migrate to many measured states and many different true states can migrate 
to the same measured state, we cannot unfold detector effects on an event-by-event basis. Rather, unfolding 
methods typically work with binned distributions. 

For the recoil library method, we chose to use a Bayesian unfolding approach [18| . This approach suits 
our needs because it is intuitive, simple to implement, and easy to extend to the multidimensional case. The 
Bayesian technique uses conditional probabilities to determine the probability that a given measured state 
corresponds to a particular true state. 

Consider a distribution of initial states E, {i = 1,2, ...,7V/}, given by P(E) (the probability of events 
with initial state ij) and a distribution of final states Fj, {j — 1, 2, Np}, given by P{Fj) (the probability 
of events with final state Fj), given the measured distribution P(Fj), and the probability for each initial 
state to migrate to each final state P(Fj\Ii), we can determine the distribution of initial states P(E) using 

N F 

P{I l ) = Y J P{h\F ] )P(F J ). (3) 

3=1 

Using simulations, we can calculate P(Ii\Fj) from P(Fj\E), the likelihood of a true state fluctuating to 
a measured state, using Bayes's theorem, which is 

= (4) 

For our particular example, with Ni initial states and Np final states, Bayes's theorem gives us 
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We can interpret this equation as follows: the probability that a given final state Fj comes from a 
particular initial state Ii is proportional to the probability density of state Ii multiplied by the probability 
that Ii migrates to Fj . The denominator normalizes the distribution. 

Our Bayesian method requires us to make assumptions regarding the distribution of initial states, P{h)- 
Although we only use P(Ii) to calculate the weights used for the measured data, the quality of the unfolding 
could depend on P(I%). To minimize this effect, the method is applied iteratively, starting with a reasonable 
prior for the distribution with Po(Ii), and with each successive iteration using the previous iteration's 
unfolded distribution as the new input. As a cross-check, the method is applied with several different initial 
Po(Ii) distributions. The iteration procedure is: 

1. Choose an initial seed input distribution for Po(7j). 

2. Using Po(Ii) and P(Fj\Ii), compute the weights P(Ii\Fj), as derived using the Bayesian equation 
shown in Eq. [5] 

3. Using these weights, recalculate the unfolded distribution P\(Ii) from the relationship P\{Ii) — 
Ef=i Po{Fj)P{Ii\Fj) described in Eq.[3] 

4. Repeat the above steps with P\{Ii) as the starting distribution. 

5. Iterate until the unfolded P(Ii) converges. 

4-3.2. Unfolding the recoil distribution 

For our application, the recoil vector is described by the coordinates (\ut\, '0'), where \ut\ is the mag- 
nitude of the measured recoil transverse momentum, and ijr is the opening angle between the measured 
recoil and the true boson direction in the transverse plane. These recoil vectors are stored in an array of 
two-dimensional recoil histograms (binned in \ut\ and ?/>*). Each histogram corresponds to a discrete bin in 
true \pt\ with bins of 0.25 GeV for small \p^\ (\p^\ < 50 GeV) and larger bins at larger |p T z |. 

In the implementation of Eq. the initial state I is specified by V^j (l^r|)fc] and the final state 

F is given by Vn; where (|p T z |)* is the true Z boson transverse momentum, (|p T z |) s is the 

smeared Z boson transverse momentum, and ip s is the opening angle between the measured recoil and the 
smeared Z boson direction in the transverse plane. 

We start with an initial seed distribution that is flat in ip*, and \ut\- We find that it takes fewer 

than 10 iterations for the unfolding method to converge. Figure Q] shows the convergence of the W boson 
mass and width obtained from the Mi distribution, as a function of iteration number in fast MC studies. 
The final value achieved agrees well with the input value. The systematic uncertainty on the W boson mass 
and width due to the unfolding procedure is discussed further in Section [5] 

Figure [5] shows an example distribution of the probabilities that a Z boson with a reconstructed \p^f \ of 

7 GeV and a recoil \ut\ of 3.5 GeV corresponds to different true \p?\ values. These probabilities are used to 

weight the given recoil as we store it in the recoil histograms corresponding to the true |p T z | . Figures ED [10] 
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show various recoil observables plotted versus the true \p^ | , obtained from the truth information of these 
MC samples, compared with the same observables plotted versus the reconstructed \Pt\, before and after 
the unfolding is applied. The unfolding corrects for average effects of \vt\ smearing on both the means and 
the RMS values of these recoil observables. 

5. Uncertainties Particular to the Recoil Library Method 

To perform high statistics tests of the efficacy of the recoil library, we study the mass and width values 
obtained by comparing |^ T | and My distributions obtained from fast MC W boson samples created 

using the parameterized recoil model with templates generated from W boson samples created using the 
recoil library method. The recoil libraries are generated from Z — > ee events created with the parameterized 
recoil method. By varying parameters in the simulation used to generate the W boson samples while leaving 
the templates unchanged, we measure the biases and statistical and systematic uncertainties on the recoil 
library method for pp collisions at ^fs = 1.96 TeV. The corresponding uncertainties for W boson mass and 
width measurements at the LHC remain to be evaluated, but are not expected to be large. 

5.1. Statistical power of the Z recoil sample 

There are significant statistical uncertainties since wc obtain the recoil system for modeling the W — > ev 
events from the limited sample of Z boson events. In 1 fb" 1 of data, after the selection cuts, we expect 
approximately 18,000 Z — > ee events with both electrons in the central calorimeter, whereas in the same 
data we expect approximately 500,000 W — > ev events with the electron in the central calorimeter. For the 
recoil library method, we choose recoil vectors from the same set of 18,000 Z — > ee events to make W boson 
templates. Our method is thus limited by the size of the Z recoil sample and any statistical fluctuations it 
contains. If we are to rely on this method as an input to a precision measurement, we need to determine 
the extent to which the statistical limitations of the Z — ► ee sample propagate to an uncertainty on the 
measured W boson mass and width. 

We assess the statistical uncertainties of the recoil method using an ensemble of 100 fast MC simulations 
resembling the statistical situation we expect in real data. We generate W and Z boson samples correspond- 
ing to 1 fb _1 of data using the parameterized recoil method. For each set of W and Z boson samples, we 
use the Z boson events to create a recoil library and then use the library to create templates for the recoil 
in the simulated W boson events. These templates are then used to extract the W boson mass and width. 
The statistical power is measured using the spread of extracted masses and widths from these ensembles. 

Figure QT] shows the measured W boson masses and widths from 100 ensembles using the Mr distribution. 
The mean fit value is in good agreement with the input value, showing that the recoil library can accurately 
model the parameterized recoil method. We test that the recoil library can model the full MC events in 
Section 6. 
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The statistical uncertainty on the mass measurement due to the recoil library method is found to be 
5 MeV from the Mr spectrum, 8 MeV for the \p^\ spectrum, and 17 MeV for the \$ T \ spectrum. These 
agree with the statistical uncertainties on the parameterized recoil method, which are found to be 6 MeV 
for the Mt fit, 7 MeV for the \p^\ fit, and 19 MeV for the |^? T | fit. The statistical uncertainty on the 
width measurement due to the recoil library method is 40 MeV using the Mt spectrum and agrees with the 
statistical uncertainty of 42 MeV using the parameterized recoil method. 

Both the parameterized recoil and the recoil library methods use the Z boson sample to model the 
recoil. One might naively expect that the additional information contained in the functional form used in 
the parameterized method would give it increased statistical power for the same-sized sample. However, 
we do not observe a loss of statistical power since the uncertainties from the two methods are very similar 
to each other. We have explored the reason for this by using a simplified detector model of W and Z 
boson events with and without recoil energy resolution effects added, and comparing the p^-imbalance (the 
difference between \pr\ and the projection of the recoil ut along the boson direction) distribution for the 
parameterized and library methods. Due to the similar transverse momentum distributions of the W and 
Z bosons, we find that the means of the p^-imbalance distribution agree with each other within statistical 
uncertainty. We also find that without recoil energy resolution effects, there is a clear but small, 0(100) 
MeV, increase in the RMS of the pT-inibalance distributions for the recoil library method, but with the 
detector resolution effects added, the RMS of the p^-imbalance distribution increases to over 2 GeV and 
masks any difference stemming from the difference between the parameterized recoil method and the recoil 
library method. 

5.2. Systematic uncertainties 

We mentioned in Section [4] that several effects could potentially bias the recoil library method. These 
include unmerged FSR photons, acceptance differences between Z and W boson events, residual efficiency- 
related correlations between the electron and the recoil system, underlying energy corrections beneath the 
electron window, and the unfolding process. The closure tests using fast MC described in Section IQl show 
the overall bias from this method to be smaller than the statistical power of the tests. Nonetheless, we want 
to make sure that this small final bias is not due to the cancellation of larger individual biases and therefore 
examine each effect independently. 

5.2.1. Unmerged FSR photons 

We measure the residual bias due to FSR photons by fitting two sets of fast MC simulations, one made 
from an unfolded, high statistics recoil file with all FSR photons included and one made from an equivalent 
recoil file with no FSR photons. We find that the mass shift between these two samples is —1 MeV for the 
M T fit, -2 MeV for the |p T e | fit, and 2 MeV for the \$ T \ fit, and that the width shift is less than 1 MeV. 
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5.2.2. Differences in geometric acceptance 

For W candidates, we only require the electron to be in the central calorimeter, while for Z candidates 
used to create the library, we require both electrons to be in the central calorimeter. To test the bias due 
to this effect, we generate two recoil files. For one recoil file we restrict both electrons to the central region 
of the detector, as we would in data. For the other recoil file, we restrict only one electron to the central 
calorimeter and allow the other electron to be anywhere, as with the neutrino from the W boson decay. We 
make templates from the two recoil files and find that the differences in both measured mass and measured 
width are smaller than the 2 MeV statistical uncertainty of this study. 

5.2.3. Efficiency related biases 

When we generate unfolded recoil files, we weight the events by the reciprocals of the ujj and SEt 
efficiencies, as described in Section [4.2.31 To check if this approach introduces any biases, we perform three 
tests. For one, we only apply the ity efficiency. In the second test, we only apply the SEt efficiency, and in 
the final test we apply both efficiencies. The maximum bias in the fitted mass or width over all three tests 
is used as the systematic uncertainty. The final uncertainty attributed to the efficiency corrections on the 
W boson mass is 7 MeV for the M T fit, 7 MeV for the |p T e | fit, and 8 MeV for the T \ fit. The uncertainty 
of the W boson width is found to be 7 MeV. 

5.2.4- Uncertainty in Aufj 

In Section 14.21 we observed that by removing the electrons from the Z — > ee recoil file, we also remove 
the recoil energy that underlies the electron cones. We correct for this effect by adding back the average 
energy, Aujj, expected beneath the electrons. We then subtract the electron energy that leaks outside of the 
electron cone that is incorrectly attributed to the recoil energy. 

We assess the systematic uncertainty due to these corrections as follows. Z boson recoil files are made 
for three cases: (1) no energy corrections, (2) a constant energy correction for underlying hadronic energy 
beneath the electron cone and constant correction for the electron energy leakage, (3) the parameterized 
energy correction for underlying hadronic energy beneath the electron cone and constant correction for the 
electron energy leakage. 

We then generate three sets of templates from each of these recoil files and measure the shift in fitted 
W boson mass and width between these three template sets. The W boson mass shifts by 2 MeV for the 
M T fit, 4 MeV for the |p T e | fit, and 1 MeV for the \$ T \ fit, with a 7 MeV shift for the width. We assign the 
magnitude of these maximum shifts as the uncertainty on this method due to the Aujj correction. 

5.2.5. Uncertainties due to implementation of unfolding 

The specific choices made in implementing the unfolding could introduce biases to the final measurement. 

Our results may depend on our choice of initial distributions in (\Pt |)*> ip> an d \ur\- They could also depend 
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on the number of iterations of the unfolding procedure we apply to the recoil library. 

We find that starting with flat initial distributions in (Ipj?!) 4 , ip, and \ut\, 10 iterations are sufficient 
to attain convergence. We generate the unfolded recoil files using 8, 10, and 12 iterations of the unfolding 
method and find that the changes in measured mass and width extracted from Mt, \Pt\i an d $t\ n ^ s are 
within 3 MeV statistical uncertainty of the fast MC study. In addition to unfolding the recoil file using 
a flat initial distribution for the recoil spectrum, we also try several smoothly varying sinusoidal initial 
distributions, and find that the variation in the final unfolded recoil file is negligible. 

5.3. Total systematic uncertainties due to the recoil system simulation 

Table [T] shows the estimated systematic uncertainties due to the recoil system simulation for 1 fb _1 
of fast MC data. The overall systematic uncertainties, obtained by adding the individual uncertainties in 
quadrature, are found to be 9 MeV using the Mt fit, 12 MeV using the \Pt\ fit, and 19 MeV using the \$ T \ 
fit for the W boson mass, and 41 MeV using the Mt fit for the W boson width. 

6. Full MC closure of W boson mass and width 

We test both the recoil library method and the parameterized recoil method using a detailed MC sample 
produced using a GEANT-based full detector model for W and Z boson production. The full MC Z boson 
sample is equivalent to 6.0 fb _1 and the W boson sample is equivalent to 2.5 fb _1 . In this case, the full MC 
Z boson samples are used to create the recoil library. Templates are then created from W boson samples 
made using the recoil library, and these are used to extract the W boson mass and width. The extracted 
values for the W boson mass and width are compared to the input values (closure test). 

Before fitting for the mass and width of the full MC sample, we test the accuracy of the model by 
comparing various full MC distributions to the fast MC model for an input value of the W boson mass of 
80.450 GeV. Good agreement between full MC and fast MC using the recoil library method is observed. 
Figure [T^] shows comparisons between W — > ev full MC and fast MC using the recoil library method for 
the Mr, \Pt\, an< ^ Wt\ distributions. The \ 2 between full and fast MC simulations are also given and 
are reasonable. The systematic uncertainties on the electron model, dominated by the uncertainty on the 
electron energy scale, are found to be 15 MeV for the Mt and f^s, anc ^ ^ MeV for the \Pt\ fit for the 
W boson mass, and 15 MeV for the W boson width. Systematic uncertainties on the hadronic model are 
taken from Section V. Since here we use the equivalent of 6.0 fb _1 of full MC Z — > ee recoils, we estimate 
the overall uncertainty due to the recoil system simulation by scaling the uncertainty due to recoil statistics 
found in Section V by a factor of l/\/6, leaving other estimated systematic uncertainties the same. The 
systematic uncertainty due to the recoil statistics is 2 MeV using the Mt fit, 3 MeV using the \pt \ fit, 
and 7 MeV using the |^ T | fit for the W boson mass and 16 MeV using the Mt fit for the W boson width, 
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which agrees with the corresponding systematic uncertainty in the parameterized recoil model. Taking the 
systematic uncertainties estimated in Section V, added in quadrature with these statistical uncertainties, 
we find the total uncertainty to be 22 MeV for the M T fit, 24 MeV for the |p T e | fit and 26 MeV for the |# T | 
fit for the W boson mass, and 36 MeV for the W boson width. 

The results of the full MC measurements agree with the full MC input W boson mass and width values 
within the uncertainties, as shown in Tabled] 

7. Conclusion 

We have outlined a method to model the hadronic recoil system in W — > tv events using recoils extracted 
directly from a Z — > 11 data library. We applied this methodology to a realistic full MC simulation of the 
DO detector. The W boson mass and width fits to these MC events are in good agreement with the 
input parameters, within statistical uncertainties. They also agree with the values extracted using a more 
traditional parameterized recoil model. Comparisons of simulated distributions using the recoil library 
method with MC give good \ 2 agreement over a full range of data observables. 

This method is limited by the statistical power of the Z boson recoil sample, as is the parameterized 
recoil model. In addition to systematic effects from the limited statistical power of the Z boson sample, there 
are several systematic effects due to the implementation of the selection efficiencies, geometric acceptance, 
the unfolding method, and FSR. The uncertainty due to these effects is found to be O(10) MeV. 

The method presented in this paper has many advantages. It accurately describes the highly complicated 
hadronic response and resolution for W boson recoils in a given calorimeter. It includes complex correlations 
between the hard and soft components of the recoil and scales the recoil appropriately with luminosity. It 
requires fewer assumptions, no first-principles description of the recoil system, and no adjustable parameters. 
At hadron collider experiments at the Run II Tevatron and the LHC, this approach to modeling the recoil 
system is complementary to the traditional parametric approach. 
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■ EM 




Figure 2: A typical W — > eu candidate as recorded by the DO detector. The two horizontal axes correspond to azimuthal angle 
and pseudorapidity, and the vertical axis is the transverse energy deposited at that location in the calorimeter. The energy 
associated with the electron and the 1ft T that corresponds to the neutrino are indicated. All other energies contribute to the 
measured recoil. The longitudinal component of the neutrino momentum is not determined, so it is displayed arbitrarily at 
Tj = 0. 
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Figure 3: Examples of the distribution of the component of the measured recoil parallel (u^ z ) and perpendicular (u±z) to the 
Z boson direction for three different bins in true \pj?\ (centered at (a) 0.4, (b) 10, and (c) 29 GeV). Each dot represents «y 
from a single event in the library. 
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Figure 4: Estimated (a) W boson mass and (b) W boson width in fast MC using the My distribution, as a function of number 
of iterations used in the unfolding. The red line indicates the input values of W boson mass and width in the fast MC. The 
default number of iterations used is 10. 
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Figure 5: The distribution of the probabilities that a reconstructed \p^\ °f 7 GeV with corresponding \ut\ of 3.5 GeV comes 
\Pt ' ' ' 
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Figure 6: Mean recoil \ut\ versus true \p^\ (black filled points) and mean recoil In^l versus the estimate of the true PjF using 
the two electrons (red open boxes) when using (a) the two smeared electrons directly and (b) the unfolded map. 
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Figure 7: Mean projection of the recoil along the Z boson direction (< u\\z >) versus true \p^\ (black filled points) and mean 
projection of the recoil along the boson direction versus the estimate of the true pj? using the two electrons (red open boxes) 
when using (a) the two smeared electrons directly and (b) the unfolded map. 
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Figure 8: RMS of the opening angle between ut and p^? versus true p^ (black filled points) and RMS of the opening angle 
between the recoil and the boson versus the estimate of the true \p^ | using the two electrons (red open boxes) when using (a) 
the two smeared electrons directly and (b) the unfolded map. 
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Figure 9: RMS of the recoil \ut\ versus true (black filled points) and RMS of the recoil \ut\ versus the estimate of the 

true \Vt\ usm g the two electrons (red open boxes) when using (a) the two smeared electrons directly and (b) the unfolded 
map. 




Figure 10: Opening angle between and true p^? (solid line) and opening angle between uj- and the estimated direction of 
true pj? (points) when using (a) the two smeared electrons directly and (b) the unfolded map for Z boson events with a true 
|p/| of 4.0 to 4.25 GeV. 
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Figure 11: (a) W boson mass and (b) width measured in 100 ensemble tests for each template generated from a recoil file. The 
dash line is a fit using a Gaussian function. All ensembles were generated with an input W boson mass of 80.419 GeV and 
width of 2.039 GeV. The fitted gaussian function for the mass has a mean value of 80.420±0.001 GeV and RMS of 0.005±0.001 
GeV. The values for the width are 2.040 ± 0.001 GeV (mean) and 0.040 ± 0.003 (RMS) GeV. 




Figure 12: Comparison plots between full MC (points) and fast MC produced using the recoil library (lines) for the W boson 
(a) Mx, (b) (c) $ t \, and (d) My (log scale) distributions. Also shown are the \ values defined as the difference between 

full MC and fast MC yields divided by the statistical uncertainty on the full MC yield. Different ranges and bin sizes are used 
for (a) and (d). 
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Table 1: Total systematic uncertainties on the W boson mass and width from the recoil library method, for 1 fb 1 of Z boson 
data. 



Source 


AM W (M T ) 


AM w (\p T e \) 


AM W (\$ T \) 


AT W (M T ) 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


Recoil statistics 


5 


8 


17 


40 


FSR photons 


1 


2 


2 


1 


Efficiency related bias 


7 


7 


8 


7 


A«j 


2 


4 


1 


7 


Unfolding 


3 


3 


3 


3 


Systematic total 


9 


12 


19 


41 



Table 2: Final result of the full MC closure fits for the W boson mass and width using the recoil library method. The full MC 
samples used here are equivalent to 2.5 fb -1 of W boson data and 6.0 fb -1 of Z boson data. For the fitted W boson mass and 
width, the first uncertainty is statistical, the second is the systematic on the electron simulation, the third is the systematic 
on the recoil system simulation due to Z boson statistics, and the fourth is other systematics on the recoil system simulation. 
AMyy represents the difference between the measured W boson mass and the input value of 80.450 GeV, and ATyy represents 
the difference between the measured W boson width and the input value of 2.071 GeV. 

Observable AM W (MeV) 

M T 6 ± 15 ± 15 ± 2 ± 7 

|p T e | 5±19±12±3±8 
\$ T \ ± 19 ± 15 ± 7 ± 8 

AT W (MeV) 
M T -5 ± 27 ± 15 ± 16 ± 10 
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